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A liquid droplet is fragmented by a sudden pressurized-gas blow, and the resulting droplets, 
adhered to the window of a flatbed scanner, are counted and sized by computerized means. The 
use of a scanner plus image recognition software enables us to automatically count and size up 
to tens of thousands of tiny droplets with a smallest detectable volume of approximately 0.02 nl. 
Upon varying the gas pressure, a critical value is found where the size-distribution becomes a pure 
power-law, a fact that is indicative of a phase transition. Away from this transition, the resulting 
size distributions are well described by Fisher's model at coexistence. It is found that the sign 
of the surface correction term changes sign, and the apparent power-law exponent r has a steep 
minimum, at criticality, as previously reported in Nuclear Multifragmentation studies [l], Q]. We 
argue that the observed transition is not percolative, and introduce the concept of dominance in 
order to characterize it. The dominance probability is found to go to zero sharply at the transition. 
Simple arguments suggest that the correlation length exponent is v — 1/2. The sizes of the largest 
and average fragments, on the other hand, do not go to zero but behave in a way that appears to 
be consistent with recent predictions of Ashurst and Holian [3, Si- 
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I. INTRODUCTION 

Fragmentation processes are important in several dif- 
ferent fields, including powder technology & H[ , soil 
physics 0, [§, p , [Ti3 | nuclear physics [ill Ha l , a stro- 
physics (IS. Il4| genetics [IH flr| . sprays Il7l. Il8l [l9j. 
fuel combustion US [SI, brittle fracture [22l. |23|. |24| and 
geology [H, [26J, just to name a few. 
Because of the nature of the fragmentation process, 
only the final distribution of fragment sizes is usually 
observable, thus becoming the main quantity of inter- 
est. Size-distributions have been characterized by dif- 
ferent empirical and theoretical statistical laws in the 
past. The lognormal distribution can be expected when- 
ever a breakdown process can be described as a multi- 
plicative cascade, and has been a ppli ed in many frag- 
mentation examples [2l|, [H, l27l. l28l. [29| . Simmons' 
law [IE HH HIJ appears to hold universally for sprays, 
and states that the mass distribution is Gaussian when 
written in terms of the square root of drop diameters. 
Villermaux [T3, [H, HH discusses the breakup process of 
ligaments in liquid jets, and concludes that each liga- 
ment gives rise to a distribution of droplet sizes that 
is essentially a Gamma distribution. The Weibull dis- 
tribution can be derived from first-pri ncip les for the 
fragmentation of solids by cracking [34 , [35j |, and reads 
%N(> m) = e -( m / A )\ where %N(> m) is the frac- 
tion of droplets with mass larger than m. Brown and 
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Wohletz (34| remark that the Weibull distribution may 
look similar to the lognormal, and suggest that the appar- 
ent goodness of fits using lognormal could be fortuitous 
in some cases. Other authors have used the (empirical) 
Rosin-Rammler form [E |34] . HE], which applies to the 
accumulated mass density in clusters of diameter larger 
than D, and states that %M(> D) = e -(D/D m )^ 
In contrast with these approaches, which postulate expo- 
nential, or else rapidly decaying size distributions, Odd- 
ershede et al [37j have claimed that the size-distribution 
of fractured objects is an inverse power-law, and fur- 
thermore sug gest ed that a form of self-organized criti- 
cality (SOC) [331 might be responsible for this. Power- 
law distributions are often found right at second order 
phase transitions, in which case one or more parameters 
must take precisely defined values. However there are 
systems whose dynamics converges to a stationary state 
that is critical, a phenomenon termed SOC [38| and for 
which no parameter tuning is needed. Later investiga- 
tions [H, EE EH EE EE Ell report "composite" power- 
laws, or transitions between lognormal and power-laws 
as the fragmentation energy is increased [5l|, [2E Ef| . 
Recent works [H, EE [43, EE EU have reported that 
fragmentation processes are only critical (and the en- 
suing size-distribution a power-law) at a precisely de- 
fined value of the impact energy. Within this pic- 
ture, the fragmentation phenomenon goes through a 
nonequilibrium phase transition at some well-defined 
value of the impact energy. Evidences of criticality at 
a well defined energy where also found in nuclear mul- 
tifragmentation (NMF) [ll], [l2j experiments, and have 
been interpreted in the context of a percolative tran- 
sition [5 01, or else in relation with liquid-gas coexis- 
tence [2|. l5ll. [52l . [53l [541 . [55l HE HE HH in nuclear matter. 
In this work, water-glycerine drops are broken by a 
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controlled-pressure gas blow, and the resulting droplets 
are counted and sized by computerized methods that in- 
volve scanning and image-processing. The only control 
parameter, the gas pressure, varies between 11 and 100 
psi, and this in turn determines the amount of energy 
available for the breakdown event. Various statistical 
measures characterizing the resulting size-distributions 
are obtained and analyzed, for a minimum of ten and 
a maximum of fifty breakdown events at each pressure. 
Our main finding is the existence of a phase transition 
around P c w 17psi, where the size-distribution becomes 
a pure power-law with an exponent T m i n ~ 1 within a 
limited size-range. Away from P c , the size distribution is 
well described by the product of a power-law and an ex- 
ponentially decaying term (Eq. (j4|)), as given by Fisher's 
model at coexistence [59| . 

The rest of this article is organized as follows. In Sec- 
tion [IT] a brief discussion is done of some theoretical 
approaches which have been used to characterize size- 
distributions in the past. Section IIHI describes our ex- 
perimental apparatus, as well as the data processing 
techniques we use to obtain statistics on size distribu- 
tions. Our results are presented in Section IIVI In Sec- 
tion IIV Al the resulting size distributions are character- 
ized by means of a simplified Fisher expression (Eq. Q), 
essentially a power law with an exponential cutoff. The 
total surface is discussed in Section IIV Bl Section IIV CI 
discusses whether the observed transition can be char- 
acterized as percolative. The mass of the largest frag- 
ment (Section IIV C II ) and moment correlations (Sec- 
tion IIV C 2p do not support such interpretation. Sec- 
tion IIVDI discusses the behavior of the largest and av- 
erage mass fragment in the context of recent theoretical 
elaborations 0, 3] . In Section IIV El the concept of domi- 
nance is introduced and applied to our data. Finally, our 
results are discussed in Section IVl 



reported to give results consistent with r « 2.3 [ill, |50| 
as well, while brittle fragmentation experiments seem to 
produce different values of r. 

NMF experiments are sometimes more generally dis- 
cussed [53l. l6lH62l 63l l63.l65l.l6r| in the context of Fisher's 
Droplet Model 59( , which states that the number of clus- 



ters (droplets, nuclei) of mass x is given by 



N(x) = N Q x- T e 



B{p,T)x-C(p,T)x° 



(1) 



where r and a are independent critical exponents. Co- 
efficients B(p, T) and C(p, T) depend on surface tension, 
temperature, and chemical potentials in the "liquid" and 
"gas" phases: 



f ' 



and 



(2) 



(3) 



Here Cq is a constant that depends on surface tension, 
and e = (T c — T)/T c measures the departure from the 
critical temperature. 

The term B(p,T)x measures the bulk contribution to 
the free energy of a droplet of size x. In a (p, T) phase 
diagram, B < if the "Gas" phase is thermodynamically 
stable, i.e. its chemical potential is lower than that of the 
"cluster", or "liquid" phase. 

The term Cx a measures the surface energy of a droplet 
of size x, so that a is the critical exponent that relates 
the dimension of the surface to that of the volume of a 
droplet. For compact droplets in 3d one has a = 2/3. 
For the liquid-gas transition a takes values close to 2/3. 
On the coexistence line, /i/ = /i g and thus B = 0, so 



N(x) = N x- T e; 



-Cx° 



(4) 



II. SIZE DISTRIBUTIONS 

In recent years, several brittle fragmentation [22l . [23l . 
[49j and NMF [111, H3| experiments have been interpreted 
in the context of Percolation Theory (6(j. A percola- 
tion transition is claimed to exist between a "connected" 
phase and "disconnected" phase, as the impact energy is 
increased. In the connected phase, a large cluster (drop, 
fragment) exists after the fragmentation event, that con- 
tains a significant fraction of the total mass. In the dis- 
connected phase no such large cluster exists. The natural 
order parameter for a percolative transition of this type 
is the average size of the larg est fragment. Measurements 
of this order parameter [23l. [50l| have shown that it goes 
to nearly zero at and above the transition, in a way that 
is consistent with the proposed interpretation. 
At the percolation critical point, the number distribu- 
tion N{x) of clusters of mass x becomes a power-law 
N(x) ~ x~ T , where the exponent r takes values close 
to 2.3 in three dimensions. NMF experiments have been 



Furthermore, at the critical point one has C = 0, and 
the size-distribution becomes a pure power-law with ex- 
ponent r. 



III. EXPERIMENTAL SETUP AND DATA 
ACQUISITION 

A. Setup 

Our experimental setup is schematically shown in 
Fig. [TJ An open, straight glass tube (T) of 6 cm length 
and 7 mm inner diameter, has a solid coaxial glass rod 
(1mm diameter) attached to its inner surface, with the 
rod's tip protruding 5 mm downwards from the lower end 
of the tube. From this tip hangs a small (11.5 fj,l [7l| ) 
liquid drop (D) of an opaque solution obtained by mix- 
ing equal quantities of glycerol and water, plus a small 
amount of water-soluble black dye. The upper end of 
this glass tube is connected to a pressurized chamber (C) 
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via a 40cm length plastic (HI) hose, and isolated from it 
by a solenoid valve (SV). The tube is fixed so that the 
droplet hangs 15 cm above the glass window of a flatbed 
scanner. 

A nitrogen cylinder is connected to the chamber through 




FIG. 1: Experimental Setup 

a plastic hose (H2), and a pressure regulator (R), which is 
provided with a manometer (P). The chamber is formed 
by a 10 cm length, 3/8 inch inner diameter plastic hose 
and is limited by two bypass valves, one manual (MV), 
the other electromagnetic (SV). 

For each droplet breakdown experiment, the nitrogen 
pressure was selected with regulator R, and the man- 
ual valve was opened until the pressure in chamber C 
equilibrated, after which it was closed again. Then, 
the electromagnetic valve was opened by means of an 
electric switch, releasing the gas contained in chamber 
C, which then flowed through tube T, fragmenting the 
droplet and projecting its fragments downwards. The re- 
sulting droplets, adhered to the scanner's glass window, 
were scanned and subsequently counted and sized with 
the help of specialized computer software. 
Before each breakdown experiment, the scanner glass was 
treated with a hydrophobic solution containing Teflon. 
Because of this treatment of the glass surface, the con- 
tact angle of the liquid onto the glass was close to 90 
degrees, so that the droplets turned out to be very sim- 
ilar in shape to semi-spheres. After the fragmentation, 
a thin spacing was inserted between the scanner lid and 
its glass window, and a high-resolution (600 dpi) scan 
was performed. Several (from 10 to 50) fragmentation 
experiments were carried out varying the jet pressure P, 
from 11 psi to 100 psi. The number of droplets obtained 
varied, from a few ones at 11 psi to tens of thousands at 
100 psi. 



B. Data Processing 

After each fragmentation, a scanning operation was 
performed. All scanned images were obtained at 600 dpi- 
grayscale, and stored using a lossy compressed graphic 



format, for later processing. Because of the non- 
negligible height of the droplets, most images contained 
faint shadows, which were removed by image threshold- 
ing 72j. Subsequent image processing and droplet iden- 
tification, count ing , and sizing, was performed with the 
help of ImageJ [73j 

The output from ImageJ consisted essentially of a set 
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FIG. 2: Data used for fitting Eq. (JS| (squares) and the result- 
ing fit (solid line). 

of droplet areas, in pixels. In order to estimate droplet 
volumes V% from the observed areas f2j, we did the fol- 
lowing. If gravity effects were negligible then all droplets 
would be plane sections of a sphere and their volumes Vi 
would be given by 



3/2 



(5) 



with Ai some constant. However (j^J is a good approx- 
imation only for the tiniest droplets, whose shapes are 
less affected by gravity. Larger droplets suffer a signifi- 
cant flattening and therefore Eq. ([5]) must be corrected. 
Under the effect of gravity, the height h(£Y) of a liquid 
droplet on a flat horizontal surface is proportional to vfi 
when Q is small, but saturates to a constant value hoc 
for large values of fi, in which case the fluid forms a flat 
"lake" of approximately constant height. Thus V satisfies 
(IS]) for small il, but is proportional to f2 for large il. As- 
suming, for simplicity, a stretched exponential crossover, 
we write 



1 



,-(n/n ) 



(6) 



where Ai,A2,^0j an d (3 are to be obtained from fits of 
calibration data (See Fig. [2]). The resulting values where 
Ai = 3.15 x 10~ 5 , A 2 = 8.3 x 10" 7 fi = 10 5 pixels, and 
j3 = 0.366. In order to obtain the data required for the 
fit, the projected areas were measured for several droplets 
having known volumes. Their volumes were obtained us- 
ing calibrated pipettes (for the larger ones), and glass 
capillaries (for the smaller ones). The data used for the 
fit, together with the resulting fit, are displayed in Fig. [5] 
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Our final expression for the volume/area relationship ex- 
trapolates to V~2x 10 _5 /d when f2 = 1 pixel, imply- 
ing that the smallest detectable volume is V m i n =0.02 nl. 
Eq. together with the parameter values obtained 
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FIG. 3: Average total observed volume Vt (squares) and av- 
erage number of droplets Ndrops (circles) as a function of jet 
pressure P. The initial volume of the droplet is 11.5 ill. 



of Xi = Vi/Y^j=i Vji t ne droplet volume fractions. The 
average number of droplets, as well as the average total 
volume Vt — Yli Vi, are displayed in Fig. as a function 
of jet pressure P. The fact that Vt remains almost con- 
stant and equal to the initial volume confirms that most 
of the mass resulting from each fragmentation gets caught 
by the scanner window, i.e. that there are no significant 
mass losses, even at the largest pressures considered. 



A. Size distributions 

Fig. 2] shows log- log plots of the normalized fraction of 
droplets N(x) with volume fraction x, for several values 
of the breakup pressure P. Continuous lines in Fig. [4] 
are the result of fitting Fisher's Model in the form in- 
dicated by Eq. Q, which is appropriate for the coex- 
istence region. We fixed a = 2/3 and adjusted Nq(P) 
(not shown), C(P) (Fig. [5]) and t(P) (the apparent ex- 
ponent, see Fig. [5]) independently for each pressure P. 
We checked that the qualitative features of these fits did 
not change much for er ranging from 1/3 to 2. 



from the fit, were embedded into a C program, which 
read droplet areas (as calculated by ImageJ), and used 
this expression to calculate volumes from them, addition- 
ally implementing the statistical analysis of size distribu- 
tions to be described in Section [TV] 



IV. RESULTS 




FIG. 4: N(x), the number of drops of relative mass x = V/Vo, 
for (from top to bottom) P =11, 15, 17, 25, 60, and 100 PSI. 
Solid lines are fits of Eq. Q with a = 2/3. Plots have been 
displaced downwards for clarity. 

For each breakdown event, a set of droplet volumes V$, 
i = 1,2, ... ,11 was obtained as described in Section Iff! Bl 
and the moments = x\ were calculated in terms 
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FIG. 5: Amplitude C(P) (circles) of surface correction term 
in Eq. (|4}, for a = 2/3. The solid line is an empirical fit (see 
text). 

The "surface" correction term C(P) (Fig. [5]) changes sign 
around Pq = 17 PSI, implying that the size distribution is 
similar to a pure power-law there, and giving the first hint 
for the existence of a phase transition. The resulting cur- 
vature change that is apparent in Fig.|¥]has been reported 
and discussed in NMF experiments as well [H, [56| . The 
solid line in Fig. [5] is a fit of the form 

<w - ,, [(s)'- 1 ] <7) 

and results in _P = 15.5 and \i = —1.5. This suggests 
that 

C~(P-P ) (8) 
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FIG. 6: Fits of Fisher model (Eq. (gj) to the droplet den- 
sity N(x) (Fig. U), produce an apparent exponent r (circles) 
showing a deep minimum T m i„ ~ 1 around 18 PSI. 



system of drops (i.e. neglecting kinetic energy) and is 
thus an interesting physical quantity in the breakdown 
experiment. We define the relative excess surface as 

E{P) = (13) 

"0 

where Qo = 1 is the initial surface. By taking aver- 
ages over breakdown events at pressure P, the results 
displayed in Fig.[7]are obtained. It can be seen in this fig- 
ure that the pressure-dependence is approximately linear 
everywhere, but there is a sharp discontinuity in its slope, 
again happening around 17 PSI. Fitting two straight lines 
intersecting at Po to the data shown in Fig. [7] (dotted 
line) one obtains the estimate Po = 17-6 PSI. 
If the excess surface is interpreted as an "energy" and 
the jet pressure as a measure of the effective "tempera- 
ture", the discontinuous slope in Fig. [7] is analogous to a 
discontinuity in the specific heat. 



near the crossing point Pq. 

The exponential term in ((4]) can be written, as is usually 
done in the scaling description of critical phenomena, as 



-Cx" 



F±((x/x )), 



(9) 



where the subscript ± means P > Po or P < Po respec- 
tively. The "cutoff" xq satisfies 

x = \C\-^. (10) 

For compact droplets we can furthermore write 

XQ = e (ii) 

in three dimensions, where £ is a correlation length (the 
linear size of the typical droplet). Considering (JSJ, (fTU)) . 
and (jlip one finds that, near Pq, 



(12) 



with v = 1/2, i.e. a classical value for a thermal critical 
exponent. Of course, this specific result (y = 1/2) de- 
pends on the validity of ([5]) , which in view of the quality 
of our data cannot be taken for granted. 
The apparent exponent r(P) (Fig. goes through a 
steep minimum T m i n ~ 1 around 18 PSI. A similar mini- 
mum in r has been taken as an indicator for the location 
of a phase transition in NMF experiments p], Q in the 
past. In fact, Richert and Wagner [TTJ] define the critical 
temperature T c as the one that minimizes t(T). We may 
thus consider this steep minimum as a further indication 
of a phase transition in the breakdown process, occur- 
ring somewhere in the range 16-18 PSI for our specific 
experimental setup. 



Excess surface 



xV 3 of the fi- 



The total normalized surface flp = JJ. 
nal set of drops in each breakdown process gives an esti- 
mation of the total amount of energy transferred to the 
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FIG. 7: Surface increase versus Pressure has a discontinuous 
derivative around 17.6 PSI. The dotted line is a linear fit, with 
the location of the kink as a fitting parameter. 



C. Is the observed transition percolative? 

1. Largest fragment 

Recently there have been suggestions that brittle frag- 
mentation [HH! EH and NMF [HI, [13] experiments can 
be rationalized in the context of a percolative [ficj phase 
transition. The validity of such interpretations rests, to 
a great extent, upon the identification of a relevant or- 
der parameter which has the properties of the infinite- 
cluster density in percolation. Within this picture, the 
connected or percolating phase must be characterized by 
the existence of macroscopic fragments, while no such 
fragments must exist in the disconnected phase (at large 
breakdown energies). The average size of the largest frag- 
ment, which goes to zero abruptly at the critical energy, 



6 



plays the role of an order parameter in brittle fragmen- 
tation [22|, HH, H^] , so it is natural to first consider the 
largest droplet size as a tentative order parameter also in 
the case of our experiments. 

The average mass Z max of the largest fragment is shown 




Jet Pressure [PSI] 

FIG. 8: Maximum droplet size Z max (circles) and average 
droplet size Z avg = M2 (squares), as a function of jet pres- 
sure P. Straight lines are of the form c\P~ 2 and C2-P _3//2 
respectively. Assuming that the system's expansion speed n 
is proportional to the jet pressure P, the observed dependen- 
cies Z max ~ P _3//2 , and Z avs — M% ~ P~ 2 are consistent 
with recent predictions of a simple energy-balance theory of 
Ashurst and Holian 0,111, 53 • 

in Fig. [5] from our droplet fragmentation experiments, ft 
turns out that Z max (P) has a smooth pressure depen- 
dence of the type Z max {P) ~ P~ 3 / 2 , without any traces 
of a singular behavior. Notice that, if there where a per- 
colative transition, the mass of the largest droplet (the 
infinite cluster) should go to zero at the critical point. 
Therefore we conclude that the behavior of Z max does 
not support the existence of a percolative transition in 
this system. 

2. Moment correlations 
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FIG. 9: Inclusive (all pressures) scatter plot of the mass Z m ax 
of the largest piece versus S2 = M2 jM\ , where the moments 
have been calculated without the contribution of the largest 
piece. Points near the "nose" of events with large values of 
S2 correspond to pressures in the range 16-18 PSI. Notice, 
however, that the existence of this feature is by itself not 
indicative of critical behavior, since it is there even for the 
trival partition described in the text (solid line). 



this, consider the breakup of a droplet of unitary volume 
into an infinite succession of droplets i = 0,1, 2,... of 
volume Vi = (1 — A) A 1 , for < A < 1. As A goes from 
almost zero to almost one, the resulting size distribution 
changes from being dominated by essentially one large 
drop, to being composed of a "mist" of tiny droplets. 
Clearly there is no phase transition in this system. 
For this simple partition, the largest mass is Z — vq — 
(1 — A). Furthermore, it is straightforward to show that 
5 2 = Z(l - Z)/{2 + Z). A plot of S 2 (Z max ) for this 
simple partition is displayed in Fig. [5] (thick line). It can 
be seen that its overall features are the same as reported 
by Campi [50j , i.e. S2 goes through a maximum and 
then decreases, as Z max decreases (A increases). Thus 
we conclude that the rounded "nose" in these plots is 
not by itself indicative of a phase transition, and cannot 
be taken as evidence of percolative behavior. 



In his discussion of NMF experiments, Campi [5(| 
notices that correlations among some moments of the 
size distribution are similar to those found in Percola- 
tion. In particular, a plot is presented of Z max versus 
5 2 = M2/M1 (moments are calculated excluding the con- 
tribution of the largest cluster) and it is argued that the 
rounded "nose" of events with large values of S2 corre- 
sponds to the percolative "critical region" . 
Fig. [5] shows a scatter plot of Z max , the largest cluster's 
mass, versus S2, for each fragmentation event in our ex- 
periments. This plot is similar in many aspects to the 
one considered by Campi (5(| , in particular with respect 
to the existence of a "nose" of events with large values of 
52- It can be easily seen, however, that this feature by 
itself is not indicative of critical behavior. In order to see 



D. Ashurst and Holian theory 

Recently Ashurst and Holian (AH) [H, H| (see also 
Astrom, Holian and Timoncn [47J), considering a d- 
dimensional fluid system that expands uniformly at rate 
77, predicted that the maximum droplet size Z max should 
behave as r\~ d l 2 while the average droplet size Z, 
should decrease 
rate 



a.s 



74] 



t?- 2 ^ 3 with 



average droplet size ^ a vg 
increasing expansion 



Since in our experiments the expansion rate 77 is not a 
control parameter, the applicability of AH results is not 
obvious. However we find (See Fig. [5]) that Z max ~ p-i-5 
and Z avg oj P~ 2 , which is consistent with AH predictions 
in 3 dimensions, if the (unknown) expansion rate rj is as- 
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sumed to be proportional to the jet pressure P. 



E. Dominance Probability 

We now set to try and identify the nature of the phase 
transition at 17 PSI, by defining a useful order parameter 
that we call dominance probability. Although, as seen al- 
ready, the largest fragment size is a smooth function of P, 
we have noticed that, at low pressures, most of the mass 
is concentrated in a few large fragments, while at larger 
pressures the mass is evenly distributed among a large 
number of small droplets. More precisely, at very low 
pressures the mass distribution appears to be composed 
of one large fragment plus many tiny droplets. Upon 
increasing the jet pressure P, often two large droplets 
appear, next events with three large fragments become 
prevalent, and so on. At large enough pressures, how- 
ever, the original droplet mass is more evenly distributed 
among a large number of droplets, without a clear size 
distinction between large and small ones. We then seek 
to define an order parameter that quantifies the property 
that the mass distribution be "dominated" by a few large 
fragments or not. 

For an arbitrarily defined "dominance factor" 7 taking 
values slightly smaller than 1 (one can take e.g. 7 = 0.9 
to exemplify the ideas), we say that there is dominance 
at level 1 if the largest fragment's mass is larger than 
7 times the total mass. If this condition is not satisfied, 
the largest fragment is removed and we check whether the 
second largest fragment's mass is larger than 7 times the 
total remaining mass (after removing the first largest). If 
this is the case, we say that there is dominance at level 
2. If not, the second largest droplet is removed and we 
proceed to do the same check with the third largest one, 
and so forth, until only the last droplet remains. 
These ideas can be formalized as follows. First all 
droplets in a given fragmentation event are ordered by 
size: x\ > x% > ■ ■ ■ > x n . Then for a given value of the 
"dominance parameter" 7, check whether the dominance 
condition 



Xk > 

k 



(14) 



is satisfied for some k < n. Stop at the smallest value 
of k that satisfies (fT4]) . In other words we will say that 
there is dominance at level kd if condition (fH)) is satisfied 
at k = kd but not at kd — 1, kd — 2, . . . , 1. If (fl"4"|) is not 
satisfied for any k < n — 1 then we say that there is no 
dominance. 

Whenever there is dominance at level kd, the first kd 
drops are substantially larger than all of the remaining 
ones. This means that the mass distribution can be sepa- 
rated in two components: one made of kd massive drops, 
the other constituted by one or more (usually many) tiny 
droplets. If on the other hand condition (fT4")) is never sat- 
isfied, it is said that there is no dominance and this means 
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FIG. 10: Probability of dominance P(k d ,P) at kd = 1 
(plusses), 2 (crosses), 3 (asterisks), 4 (empty squares), and 
8 (full squares); mass of dominant drops Md{P) (empty cir- 
cles), and total dominance probability Pd(P) (full circles) as 
a function of pressure P. 



that the distribution has no distinguishable size scale. 
In NMF, events where the excited nucleus emits a few 
small fragments composed of one or two nuclei, remain- 
ing otherwise almost unaltered, are called evaporation 
events. Those events where the original nucleus splits 
into two large fragments (plus eventually a few much 
smaller ones) are called fission events, and those for 
which three large fragments are formed are called ternary 
fission events. Within the dominance concept we have 
introduced, the NMF classification mentioned above cor- 
responds to kd = 1, 2 and 3 respectively. 
In order to apply these ideas to our experimental data, for 
each pressure P we determine the probability Po(kd,P) 
that there is dominance at level kd, by averaging over 
fragmentation events (typically fifty) at pressure P. The 
resulting dominance probabilities are displayed in Fig. 1101 
for several values of kd- It is apparent from these data 
that evaporation (kd = 1) events are predominant at 
low pressure, next fission (kd — 2) events appear as 
P is increased, followed by a predominance of ternary 
fission^ = 3) and so forth. As shown in Fig. [ITJ the 
average dominance level kd grows on approach to P = 18 
PSI, beyond which there is no longer any detectable dom- 
inance. 

The sum over k of Po(k, P) gives probability D(P) that 
there is dominance at some level, and is also shown in 
Fig. [3Z3 together with the average mass Mfj(P) of all 
dominant droplets. The fact that Mfj(P) is roughly 
equivalent to D(P) indicates that, whenever there is 
dominance, almost the total mass of the system is con- 
tained in the kd dominant drops. The total dominance 
probability D(P) goes to zero at around 18PSI, in a 
way that is consistent with all previously mentioned evi- 
dence indicating a critical point around that pressure. Of 
course the precise numerical values of the dominance ob- 
servables Po(k, P) depend on the arbitrary choice of the 
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FIG. 11: Average number of dominant drops, versus pressure. 
For pressures larger than 18PSI, no dominance was ever de- 
tected. 

dominance parameter 7, but we have checked that the 
fact that the dominance probability goes to zero around 
18 PSI holds qualitatively unaltered for 7 ranging from 
0.8 to 0.95. We also explored a different definition of 
dominance, namely replacing (fTi]) with 

x k > jXk+i- (15) 

This means that a comparison is established between 
each mass-sorted droplet and the next, stopping at the 
first level at which the ratio between them is larger than 
a given prespecified factor. The use of definition (|15p in- 
stead of (H3]) did not alter our findings qualitatively. 
We thus feel justified to conclude, based on the data pre- 
sented in this section, that there is a dominance transi- 
tion in the neighborhood of P c = 18 PSI. On the low- 
pressure side of this transition, a few large (dominant) 
drops contain most of the original mass, and there is a 
size-gap separating these large drops from the accompa- 
nying "mist" of tiny droplets. Beyond P c , on the other 
hand, one has a collection of droplets which populate 
homogeneously the whole range of sizes without a sys- 
tematic size gap. 

V. DISCUSSION 

We have presented a statistical analysis of droplet-size 
distributions obtained when a liquid drop is broken by a 
sudden gas jet blow. A flatbed scanner and image pro- 
cessing software were used to automatize the process of 
counting and measuring the resulting fragments, enabling 
us to collect statistics for breakdown processes resulting 
in up to tens of thousands of tiny droplets. The result- 
ing size-distributions N(x), displayed in Fig. |H can be 
fitted by a simplified Fisher- model expression (Eq. ((4])) 
and it was found that N(x) behaves like a pure power- 
law only at P c w 17 PSI, providing the first hint for the 
existence of a phase transition there. Further evidence 



suggesting a phase transition is given by a deep mini- 
mum in the apparent Fisher exponent r (Fig. [5]), and 
by a sharp slope discontinuity in the relation between 
pressure and surface increase (Fig. [7|) , both of them oc- 
curring at P c . We argued that it is not possible to iden- 
tify P c as a percolation critical point, as was done in 
NMF p], El and Fracture [H, H| [1| experiments, be- 
cause: a) T mm w 1 (Fig. [5]) is far from the value 2.3 
that would be expected for a percolative transition, and 
b) in our case, the mass Z max of largest fragment does 
not become zero at P c , but instead behaves as P~ 3 / 2 in 
the whole pressure range considered (See Fig. |5J). Low 
values of the apparent exponent r have been previously 
reported in NMF experiments fll H. 167} . 
Furthermore, the behavior of the largest and average 
fragments, in our experiments, was found to be consistent 
with recent predictions of Ashurst and Holian (AH) [HIH, 
if the original droplet is assumed to expand at rate r\ that 
is proportional to P (Fig. HJ). 

Analysis of moment correlations through a scatter plot 
(see Fig. |9]) reveals a behavior similar to the one already 
found by Campi [50I ] in NMF experiments, and which 
was invoked as evidencing percolative critical behavior. 
We showed, however, that the shape of the scatter plot 
displayed in Fig. [9j i.e. the existence of a maximum in 
the values of S2(Z), is in no way related to criticality, 
since it is there even for a simple, non-critical, partition- 
ing model, as discussed in Section HV CI 
We introduced the concept of "dominance" in an attempt 
to quantify the observed fact that a few large droplets 
contain most of the original mass at low pressures, while 
no dominant droplets are observed at larger pressures. 
The dominance probability, as we defined it, although is 
dependent upon an arbitrary parameter 7, shows a sharp 
fall (Fig. fTU|) to zero at a critical pressure P c w 18 PSI, 
that is roughly independent of 7. Of course the precise 
value of P c is not interesting by itself, since it will de- 
pend upon many parameters like the volume and viscos- 
ity of the original droplet, the amount of gas contained in 
the pressurized chamber, the diameter of the outlet (and 
thus the exhaust velocity), etc, whose variation we did 
not consider in this experiment. Conceptually more in- 
teresting is the fact that, as the energy available to break 
the droplet increases, a sharp transition occurs between 
two breakdown regimes, and that this transition can be 
characterized by statistical methods as discussed in Sec- 
tion [Tv] 

Future work could include the use of fast imaging tech- 
niques in order to put these results in the context of 
what is known about rupture modes in the field of 
sprays [H, [68|, [69|, [7(|. The following picture emerges 
for isolated droplets subject to a sudden gas blow. At 
low excitation intensity (low Weber number) the preva- 
lent rupture mechanism is vibrational breakup, whereby 
collective oscillations break the drop into a small num- 
ber of roughly equal-sized droplets. Upon increasing the 
breakdown energy, bag breakup appears, which is char- 
acterized by the deformation of the drop into a torus- 
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shaped rim, which subsequently disintegrates. A variant 
of this is called bag-and-stamen breakup, in which a cen- 
tral "stamen" is also formed. Next come shear breakup, 
where small droplets are continuously stripped off the 
rim of the drop, and catastrophic breakup, where strong 
surface waves disintegrate the drop violently. The limit 
Weber number between vibrational and bag breakup is 
around 10. At 20-60 shear breakup appears, and at 1000 
catastrophic breakup. 

Since our experiments were conducted on droplets hang- 
ing from a thin glass tip, the presence of which certainly 
modifies the gas flow and the ensuing breakup modes, 
we do not expect a direct correspondence with the above 
classification, although a comparison with the rupture 
modes of free droplets would certainly be instructive. 



Also interesting is to establish a comparison of our re- 
sults with recent analyses of NMF data in the context of 
phase transitions [HI, [f| • Common characteristics are: 
a sharp minimum in the apparent r exponent [j], Ip. [Till . a 
strong curvature change in a log-log plot of A^(%)[53|, [56j, 
and the fact that N(x) is, to a good appro ximation, a 
pure power law at criticality (ll. SflTI. lEa. l56|. 
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